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ABSTRACT 

Particle acceleration in relativistic shocks is not a very well understood subject. Owing 
to that difhculty, radiation spectra from relativistic shocks, such as those in GRB after- 
glows, have been often modelled by making assumptions about the underlying electron 
distribution. One such assumption is a relatively soft distribution of the particle en- 
ergy, which need not be true always, as is obvious from observations of several GRB 
afterglows. In this paper, we describe modifications to the afterglow standard model 
to accommodate energy spectra which arc 'hard'. We calculate the overall evolution 
of the synchrotron and compton flux arising from such a distribution. We also model 
two afterglows, GRB010222 and GRB020813, under this assumption and estimate the 
physical parameters. 
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1 INTRODUCTION 

Relativistic particles accelerated by shocks occupy a predominant place in astrophysical systems. These particles emit syn- 
chrotron and compton radiation, which can be observed from radio to gamma-ray bands. Gamma Ray Bursts (GRBs), their 
afterglows, Supernova Remnants (SNRs) , Active Galactic Nuclei ( AGNs) and Pulsar Wind Nebulae (PWN) are some of the 
most important and intriguing candidates which house shock accelerated electron population. 

The details of these electron populations and hence the details of the acceleration process are inferred from studying the 
emitted synchrotron and compton radiation. The accelerated particles are often found to be distributed non-thermally, as a 
power law in energy characterised by an index p : 



iV(7e) = ii'eTe ^ (7m ^ 7e <lu), 



(1) 



where N{'y)d'y is the number density of electrons in the energy interval ■yrUe C^ and (j + d 'y)meC^ ■ 

This non-thermal power law is a natural outcome of the Fermi process ( Fernij[ 194E ), a standard framework to describe 
shock acceleration. Several analytical and numerical investigation s have been made ( Achterberg et al.ll200ll : IOstrowski fc Bednara 
2002l : lElhson fc Doublj|2004l : lKeshe^l2006l : lNrshikawa et aLlbood ) especially for the Diffusive Shock Acceleration (DSA) mech- 
anism, a variant of the Fermi first order process, which is expected to operate in coUisionless shocks. Most of the theoretical 
and numerical studies produce a 'single soft' distribution of the accelerated particles, where the index p is greater than two. 
Though there are many observations supporting this prediction, a non-negligible fraction seems to differ from this. 
Observations of some AGNs , GRB Afterglow s and PWNs have revealed an underlying 'hard' (p < 2) electron distribution 
( Panaitescu fc Kumaiil2001al : IShen et al.ll2006l ). Derivation of expressions for the radiation spectrum from such a distribution 



requires a different treatment from its 'soft' counterpart. 

In this paper, we introduce a modelling platform for afterglow spectral evolution in the presence of a hard electron 
(p < 2) energy distribution. We then present the model of a few afterglows with such a hard spectrum, and derive their 
physical parameters. 
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2 HARD ELECTRON ENERGY SPECTRUM 

The distribution described in equation- 1 can be safely assumed to go to infinity if it is soft, since the role of the higher energy 
end is negligible in total number and energy content of the distribution. Hence the equations which form the basics of the 
standard afterglow modelling paradigm contain only 7^ and p. 

However, a hard electron distribution can not be extended upto infinity, and requires to be terminated with an upper 
cut off to keep the total energy from diverging. This upper cut-off, 711, which is determined by the acceleration mechanism, 
plays a crucial role in the analytical treatment of p < 2 spectra. Since electrons towards the higher energy end dominate in 
the share of the total energy content in the distribution, the upper cut-off appears explicitly in the equations describing the 
spectral parameters. The distribution beyond 711 could be a sharp drop, an exponential fall off, or a steeper (p > 2) powerlaw 



There have been previous studies to incorporate hard electron energy distributions in afterglow modelling. iBhattacharva 



2OOII (hereafter BOX) has used a 7^ which is a function of the bulk lorentz factor (F) of the shock. The dependence on F is 
parametrised by an index q. 

The time dependence of 7m is altered by the introduction of 7„. This in turn modifies the spectral evolution. Moreover, 
a n ew break frequency corresponding to 711 will appear in the spectrum. 

Dai fc Cheng|l200ll (hereafter DCOl) has followed the same approach but by constraining 7^ (in their notation, "/m) to 



be due to the termination of acceleration process by energy loss to synchrotron radiation. Their model is a special case of lBOl 
wit h q = —1/2. This upper limit 7m, in typical conditions lie at very high energies 



Panaitescu fc Kumaij l2001bl (hereafter PKOl) consider two conditions to determine the upper limit of the energy dis- 
tribution, (i) The upper limit (7M1) results when the acceleration timescale becomes larger than the timescale for radiative 
energy loss (same as DCOl), and the corresponding break frequency lies much above the observation limit, (ii) In the second 
case, the distribution terminates at an upper cut-off (7A/2). A steeper powerlaw is assumed beyond the cutoff. A constant 
fraction of the shock produced thermal energy is assumed to be contained in the electron distribution, the lower bound of the 

distribution 7™ is assumed to follow the same evolution as it does in the standard model. The evolution of 7M2 results from 

p-i 
these two conditions. In the limit, 7^/2 3> 7™, and F 3> 1, 7A/2 can be obtained analytically to be proportional to F p-'^ . 

The second assumption that 7™ follows its standard model behaviour is somewhat inappropriate in this context, since this 

behaviour corresponds to a condition where the effect of 7A/2 is ignorable. In reality, 7A/2 originates in some physical process 

which will have its own dependence on F, hence it is more appropriate to paramet rise t he evolution of 7A/2 as a function of F. 

In this paper, we continue the investigation of iBOll . The upper cutoff 7„ of iBOll is identified as an injection break 7^, 

above which the electron distribution steepens to a powerlaw with index p2 > 2. We leave room for accommodating different 

processes, by keeping the parametrisation of 7^ to be that of BOl. Our results differ from PKOl in having the evolution of 7™ 

and hence of the lightcurve, depending on the nature of the injection break. The flux decay index and the closure relations 

between the lightcurve decay slope and spectral slope also depend on the injection break, essentially the value of q, which is 

characteristic of the mechanism responsible for the upper cut-off. 



3 MODIFIED ELECTRON DISTRIBUTION AND INJECTION BREAK 

The double slope electron energy distribution with slopes pi and p2 is represented as. 



NM 



Keje''^ 7m ^ 7<: < 7ii 
K'e^e^'^ 7i ^ 7e < 00. 



(2) 



Here, K^ is the normalisation constant, which will depend on the number density of the ambient medium n{r-) and the bulk 
lorentz factor F. K'^ can be written as, 7^6 7^. 
We modify the BOl parametrisation of 7i to 



7. = C(/3r)' 



1 ^ 7i ^ 00. 



(3) 



in order to accommodate the non-relativistic regime of expansion where F ^ 1 and /3 9^ 1. Usin g the standard r esult that the 
post shock particle density and energy density are 4Fn(r) and 4F(F — l)n{r)mp(? respectively ( Sari et al.lll998l ), one derives. 



A'e = 4n(r)gp ^ 



me C^-Pi /J-JCa-Pi) 



[F - 1] F 



[l-<7(2-pi)l 



(4) 



7m = 



me i 



2-pi 



u 



(pl-1) 



1 



piCi-Pi) 



(pi-i) 



[r- 



1 l(2-Pl) 

i]¥r^ F pi-i 



(5) 



where, rup and rrie are the proton and electron rest mass respectively. The function gp = fp(pi — 1) and fp = -r — _]'?'■},_ — K- 
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3.1 New Spectral Break 



The standard afterglow model has four spectral parameters, the synchrotron peak frequency, Urn, the cooling break or the 
synchrotron cooling frequency, Uc, corresponding to the lorentz factor beyond which the electrons cool rapidly, the flux /p at 
the peak frequency {vm or fc), and the synchrotron self absorption (SSA) frequency, Va, above which the fireball is optically 
thin. The radiation spectrum emerging from a double slope electron distribution will exhibit an additional " injection break " , 
corresponding to the lorentz factor 7i. Using the standard expression for synchrotron frequency i^syn(7) for an electron lorentz 
factor 7, one obtains, 

1 + z nnieC 
where B is the post-shock magnetic field density, e is the electron charge, c is the velocity of light and z is the redshift of the 
burst. 

Above this frequency the spectral slope steepens to the value corresponding to p2 from that of pi. 



4 SPECTRUM : THE SOURCE FUNCTION METHOD 

Instead of the usual approach of writing fiux f^ en u^ , we use the synchrotron source function along with the optical depth 
to obtain the final flux. Therefore, 

./. = S^l - exp(-r,.)] (7) 

where 5*1^ is the synchrotron source function, which has the following functional form: 

2 

ly < i^p 

,5/2 (8) 

" ^ l/p 




u > u„ 



Svp is the source function at peak frequency Vp. For slow cooling (ie., Vrn < Vc), Vp = Vm and for fast cooling (ie., Vc < Vm), 
Vp = Vc- Svj, can be calculated as, Sv^, ~ fpTvp, where fp is a normalisation constant, that equals the fllux that would have 
been expected at Vp if self absorption were absent. 

The optical depth due to synchrotron process varies as j^~^" when v is less than Vp and i/~(p+*"2 otherwise. Normalising 
the optical depth to be unity at i^ = j^a, r^p, the optical depth at u — Up can be written as [^1 when Ua < Vp and 

[^1 when Va > Vp- Value of p in the latter expression is 2 for the fast cooling regime \i Vc < v < Vm- p is replaced 

by Pi (pi + 1) and p2 (p2 + 1) in the slow (fast) cooling regime below and above Vi respectively. 
For a double slope electron energy spectrum undergoing slow cooling, 

-5/3 

V < Vm 

-(pi+4)/2 



X < 



{-;;-) Vm <v < \Vc,Vi) 

/ JI\-(Pl+4)/2 / ,^ \-(pi+5)/2 ^ ^ ^ CQI 

(U^j (it) V^<Va<V <Vi (,yj 

/ Jr\-(Pl+4)/2 / ^%-(P2-l-4)/2 ^ ^ ^ 

KTt) \-) Vrn<V,<V < Uc 



For fast cooling, 

( V \-5/3 



V <Vc 

Vc <V < V„ 



3 

(10) 



(vr„\^^( u \-(Pl + 5)/2 _, ^ _, 

\'^) VI^/ Vc<Vttx <V < Vi 

[^) [it) (-) Va<Vn.<V,<U 

These expressions, along with equation-|8] are substituted in equation-!?] to obtain the final flux, which at a given time, is 
a function of the five spectral parameters {vm, Va, Vc, Vi and /p). 



To estimate these parameters, we first evaluate r(r) and r(f). For that, we use the expressions given bv iHuang et al. 



( 200(J ). after correcting for redshift, which accommodates a smooth transition from an initial ultra-relativistic to the final 



non-relativistic regime of the fireball. Time evolution of the half opening angle [Qj) depends on the lateral velocity of the 
jet in its comoving frame, which essentially is the sound velocity of the post-shock medium. The half opening angle varies 
as, -^ — -^ [if]' where Cs is the velocity of sound in the downstream medium. Cs is usually assumed to be constant 
throughout the evolution of the shock, but this is not a very accurate assumption. Initially, when the downstream plasma is 
ultra-relativistic, the thermal velocity will be c/\/3, but as the ejecta becomes non-relativistic, the velocity approaches \j — 
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where lUp is rest mass of the proton. We calculate Cs as a function of F, adopting the method followed by IChandrasekhai 
( I939I ). This gives us 



ksT 1 



(11) 



We have used equation- A3 (Appendix-I) to obtain temperature in terms of F. More details of the calculation is given in the 

Appendix. The comoving magnetic field density B is given as [STrea y 1 ^ c, where es is the fraction of thermal energy in 
the magnetic field, m is the total swept up mass, 14o is the volume of the downstream plasma in the comoving frame, which 
can be calculated as Qr^A' where Q, is the solid angle and A' is the comoving shel l thick ness. 

We calculate fp using the expression (equation-25) given bv lWiiers fc Galamal ( 1999 ). Um and Vc are calculated using the 
expression described in section IHTTl by replacing 7^ with 7™ (equation-[S]) and 7c {— &Tvmec/{aTVB^t)). Va is the frequency at 
which the synchrotro n optical depth in the comoving frame (a'j^, A', where a,j is the absorption coefficient calculated following 
the method given bv lRvbicki fc LightmarJ ( 19791 ')') equals unity. 

For various values of q, the evolution of the spectral breaks as a function of time is plotted in figure 1 and the lightcurves 
are displayed in figure 2. The difference of evolution introduced by q is apparent in these figures. 



5 DYNAMICS : LIMITING CASES 

To obtain the overall dynamics of the fireball, we adopt the method presented bv lHuang et al.l (2000|) which accomodates a 



smooth transition from the initial ultra-relativistic to the final non-relativistic phase. 

However, analytical solutions for F(r) are possible in extreme cases. The adiabatic (e = 0) ultra-relativistic regime 
(F ^ l,/3 ~ 1) is encountered most commonly in afterglow observations. At late times, {t > ^nr, the fireball becomes 
non-relativistic. This phase is same as that of the well studied supernova remnants. 



5.1 Ultra-relativistic Limit 

In this limit, the expressions for F(r) and r{t) of lHuang et al.l (2000|) can be approximated to 

yjii- s)Eo/{nc^) {parlY^^'^ (r/ro)'""^'/^ and ((4-s)(3-s)2cf£o/((l + 2)!^c^Poro))^ respectively, where p(r), the ambient 
medium mass density profile is parametrised as po (r/ro)'" ■ The expressions for spectral parameters we obtained for this phase, 
are listed below. We consider two types of ambient media, (i) a constant density around the progenitor star {n{r) — n,s = 0) 
and (ii) a stellar-wind blown stratified density profile (s = 2, with a normalisation po = 5 x 10^ j4* and ro = 10^'^ cm). 



fp(mJy) 



210.45 



1021.5 %^i±^ 

4,Gpc 




yj£iao,h2(-BA 



l-9(2-pi) 



(s = 2) 



1.87x10^17.14) ^1-1 r^/Jpi-i ^^^J^^e. 



Vm{Hz) = < 



C" 



,2-pi 



£i3o,52 



pj+qpj— 2g 
■SlPl-l) 



C 1 + Z 
3(Pl+gPl — 2q) 



(1+^) 



5.77 X 10^ (13.1)- yf ^ [^U.]^^\ £^,2 Ai'-'^^' 



trf 



(1+z) 



where y 



Va{Hz) 



MHz) = 



l-q(2-pi) 
Pl-1 ' 



bp and Xp are functions of p ( Wiiers fc Galamalll999l ) 



5.84x1013 [1^] C' 



1/2 



-1 -3/2 



7.6 X 10' 



1.3 X 10" 



(1+^) 



Cs -3/2 ^ . 



[U{l + z)]-'/^ (s = 0) 

1/2 



-3/2 .-2 0I/2 

^* '^iso,52 



2) 



£iso,52 



1.65 X 10^ ^\3.62)''y^£^i%Ai^+'''>^\]^^ 
In the slow cooling regime. 



1(1 + 9) 



TT+^ 



'^{2+q) 



-(1 + 9) 



(8 = 0) 

(s = 2) 



(s = 0) 
(s = 2) 



(12) 



(13) 



(14) 



(15) 
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^a(Hz) 



2(P1-1) 



2+Pl 



Pl+1 



(2.61 X 10^^)''+pi (8.38 X 10")-i+pi (8.2x10"^) -i+pi (0.88) "i+pi (3)pi+* 

PI +2 
2(2-2<; + pig) 2 + pi r ^ -1 I^TJ 

(4.14) *+pi (0.64) *+pi y^ ' 

(2.3 x 10-''')5^ [r(^) r(S^)]^/^"+*' [1.87 x IQ-^^gpVS^] ^ 

Pl+6 + gpi-2<; Pl+2 pi+6^opi+2n 2- 

4(pi+4) 2(pi+4) 2/(p,+4) ,, ''j^^^^'' ^-2^-3- 

^iBo52 <:S Ee n -"(Pl+I) ^ P1 + 

— 10 — 3p2^ — 3gpi +6g 
\Jd_] 4(pi+4) 

Ll + Z I 



PI 



(S = Ua > Vm) 



lQg4-3pi — 5gpi - 

2.61 X 10" (4.14) 5(Pi-i) 

3/5 



1 23/40 



(4-pj)(P2-2) 



-1 
1/5 ,(Pi-l) 



yp m. 



18-13pi-fQg + 5piq 



(pi+2/3)(p2-Pl) 
3(Pl-2)(g-l') 



where T denotes the Gamma function. 



S ISO, 5. 



40(pi-l) 



14-19pi+lQg-5pig 
J, 40(pi-l) 



(16) 



I^a(Hz) 



(3.62)^2 (3.42 X 10'''^)(pi+'') (2.4 x 10^)(pi+'») 



p2(pi+4) . 
•^130,52 ^ 



(Pl+4) 



[v/f^^¥]"+^ r(M±2) r(^ai±H) [£.]5TT5 

i;Pl-2.; l. + 2pi+4g-2piq ^ , 2g-2pi -gpi- 

FTT3 



6.16 X lO"'* 3.62 



14 



1H3 



Pl+2 



(1+z) 
3/2 



5p J2-P1 



(s = 2 Z/a > Z^77i) 



3/5 1/5 r.5 
?P ^S '--iso,52 



pi+2/3 

2 I (Pl-2)(l-<;) 6 (Pl-2)(l-g) r 
4(P1-1) /l S 



2 + 3pi 
4(P1-1) 



(1 + 2) 



7pi-2-10g + 5pii; 
20(pi-l) 



52-Pi 



5(P1-1) 



(S = 2 Va < Vm.) 



where 2/2 = Pi+s-^g+^gpi and yg = '"^'^"''"^ 

^ pi+4 ^^ PI— 1 

In the fast cooling regime, 



-1/2 



i^a(Hz) 



1.96 X 10" (pi - l)^/« [^]'^^ (1 + zf/' £lj:^,, n^' [if 

1.83 X 10- (PI - 1)^/^ [^] '"'' (1 + .f £,1^52 nj^^" 4^^ [ife] "^'' 
1.44 X lO'-' (pi - 1)'/=* AV" [U{1 + z)]-^''-" 



(17) 



1 6/5 



8.48 X 10^ (pi - 1)=^/^ \^r^ £i 



2/5 .11/5 6/5 



M^ 



8/5 



J^c < J'a < I'm (S = 0) 

I^a < J'c (S = 0) 

J'c < J'a < I'm (S = 2) 

i^a < I'c (S = 2) 



(18) 



5. -Z.i a-5 closure relations 

The a-(5 closure relations for a general value of q valid in the slow cooling phase of the ultra-relativistic approximation are 
the following: 



f [{q - 1) - 25{q + 1)] v„, <v <Uc, t< tj 
i[(3g-l)-35(q + l)] v>Uc, t<tj 

Vm < V <Vc, t > tj 
V > Vc, t > tj 



\[q-?,-25{q + l)] 
{q-l)-5{q+l) 



(19) 



In figure 3., we display the above closure relations. The q = 1 plot can be considered as a reference to the standard model, 
as it recovers the usual slopes. The dependence a has on q has to be kept in mind while inferring the value of p from the 
lightcurves. Temporal decay indices calculated for the ultra relativistic limit are listed in table 1 and lightcurve decay indices 
are listed in table 2 (slow cooling) and in table 3 (fast cooling). 

5.2 Non-relativistic Limit 

In the non-relativistic limit, at i = ^nr, the lorentz factor is ~ 1 and /3 ^ 1. The fireball by this time would have undergone a 
considerable lateral spread and the geometry may be approximated to be spherical. The solid angle SI may now be set to 47r. 



5.2.1 Dynamics 

The evolution of the radius r is calculated as, 
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r = (^(7) 



Eo 



Po 



-|l/(5-s) 



^2/(5-s) 



(20) 



where Eo is the energy in the explosion and 7 is the ratio of specific heats for the plasma. One could assume (^(7) to be 1.05 
for a constant density ambient medium and 0.65 for a stellar- wind blown medium ( Berger et al.ll2004h . 



5.2.2 Electron energy spectrum 

The thermal energy density in the shock downstream is estimated as, 



Wth 



gc'po 



/3' (r/ro)"^ 



where /3 is \^. The expressions for electron number and energy will give, respectively. 



K, 



_^ = 4po/mp (r/ro)" 

1 - IJTm 



KeTUe (2-pi) 9c^P0 n2 , , ^ 



(Pi - i)t 

5p ^' 2 

Solving eq. [23] and eq. 1221 one obtains the expressions for K^ and 7„ 

^^='5p^^(r/ro)-»/3^-(^--' 



2"" me e 



7m = 



^ ^ iii'p ce 



i/(pi-i) 



/3 PI 



(21) 

(22) 
(23) 

(24) 
(25) 



5.^2.5 Spectral Parameters 

The magnetic field energy density is assumed, as usual, to be a fraction es times the thermal energy density, ie. 



B = v/97reB/?2c2p(r) 

We calculate the four spectral breaks, Ua, i^m, i^c, and Ui and the peak flux fp from: 



(26) 



/. = '•''." '""" 0.053"'^-^> r(^) r(^Ei+2) r^:^ 



dla 



v^ = 2.8 X 10*^ -^ B 

(l + z) 

V, = 4.81 X 10^^ -^ 
i/i = 8.0 X 10^ 



iNR 



2065.7/,^_^^, 



iNR 



/3 (pl-1) 



1 + 2 



(1 + -) 



5r/3 



2fl2g 



4.72 



Pl+2 
Pl + 2/3 



K, 



3/5 



^_(3p,+2)/5 ^3/552/5 (fo^ ^^ > ^^) 



J^a = < 



(6.72 X 10"^'^)pi+-i (1.25 X lO^'')""-' (7 x 10~'')pi+'' 



(27) 

(28) 

(29) 
(30) 

(31) 



I [r(M+H) r(^El±H)]5IT4 BiT+4 [rS^] pi+4 (for V, < v^) 
where a is a numerical factor, describing the thickness of the shock in terms of r as A = r/a 



6 SYNCHROTRON SELF COMPTON EMISSION 

The contribution to the total flux from synchrotron photons which are compton scattered by the non-thermal relativistic 
electrons themselves, can be signiflcant towards higher energies. ^^^_^^^^ ^_^ 

We calculate this compton component following the method adopted bv ISari fc EsinI ( 200 ih . Following this work, the 
approximate ratio of inverse compton (IC) to synchrotron luminosities may be estimated as follows (for a uniform density 
ambient medium and the slow cooling regime). 

The IC spectrum is characterised by four break frequencies : v]^ = 27^1^^", v]P = 2'y'^v^^^ , v] = 2'yfu^^''^, u^ ~ 27^;^^^" 
and a flux normalisation fp = fp^'^ur nr 



Hard Electron Energy Distribution 7 

For !/^" ^ u^^^^ ^ i^c^"'i the energy emitted by the compton process peaks at ul'-' and that by the synchrotron process 



will peak at u^' 



IC rlC 






, syn rsyn 
'^C J Da 

7007^_7 7c,7 



(32) 



7m, 500 



" (pi-l)o.; 



7i,5 



7i,5 



" (P2-l)l.5 



7c,7 



For !/S 



^ 1/=^" $: 



compton energy peaks at vf^ and synchrotron energy peaks at j/J* 



IC i^IC 



, syn rsyn 
7007^-7 7i,T 7c,5 



(33) 



7m, 500 



where 7e,i7 



7i.7 
7e/10", 7m, 500 



(pi-l)o.5 






and {p — 1)/ = (p — I)//- In either case, the compton power 



7™/500, 71-7 -^o_ 

peaks at very high frequencies, (~ 10^^ Hz g-fQ y|jo) for si, hard electron spectrum. Hence the contribution of synchrotron self 
compton emission becomes significant only at frequencies above hard x-rays. 

As a next step, we estima te the IC flux from a numerical integration over the photon and the electron spectra. To do so 
we use the expression given bv lSari fc EsinI (|200ll ) for the inverse compton flux due to the modified electron distribution, and 
the synchrotron radiation spectrum /^■*'" generated by this electron energy spectrum. 



J.IC 

U = rcTT 



d-jNi-j) I dxf7^{x) 

'0 



(34) 



where xq ~ 0.5 

The synchrotron and the compton fluxes obtained from the above calculation are displayed in figure 4. 



7 MODELLING SHALLOW EVOLUTION 

A new parameter q is required for modelling afterglow evolution based on hard electron energy spectrum. This index 
parametrises the evolution of the upper cut-off of the electron spectrum (see equation-3). The value of q is determined 
by the acceleration process operating in the relativistic shocks. The present understanding about this from theoretical or 
numerical calculations is not exhaustive. 

The termination of the acceleration process due t o synchrotron radiation losses leads to 7i being in versely proportional to 
the square-root of the bulk lorentz factor {q = —0.5) ( Gallant fc Achterberdll999l : lLi fc Waxmanll2006l ) . However, the slowest 
post jet break decay in this case tends to 1.75 as p\ tends to its minimum possible value of 1 (in the limit 1 ^ pi ^2). This 
is noticed by DCOl also, who have tried to model GRB010222 using a hard electron energy spectrum. They have used this 
fact to rule out the presence of a hard electron energy distribution in this afterglow. None of the afterglows we model in this 
paper however display post jet break decays steeper than 1.75, which rules out the possibility of their electron distribution be 
terminated by synchrotron losses. 

q = 1 is applicable to the lower cutoff of fermi pro cess (7^ = ^^^ F), below which a pre-acceleration mechanism producing a 
flat electron spec trum may operate ( Achterbereil200ll) . The presence of such an u pper cut-off is observed in som e of the Active 
Gala ctic Nuclei ( Leahy et al.l 1 19891 : iKonopelko et al.l 120031 : IStawarz et al.l [20071) and Pulsar Wind Nebulae ( Hoshino et al. 
19921 ). Moreover, g = 1 also provides scalings that would have been obtained in the standard fireball model without references 
to 7i . Good fits could be obtained with a g of 1 for all three afterglows we study (JBhattacharva fc Resmill2004l : iMisra et al. 
2005|), however, the value of ^ we inferred from these fits are far higher than mp/m^. 

Another interesting value of q is —1.0, though any mechanism producing such an upper cut-off proportional to the inverse 
of the bulk lorentz factor is not discussed in the literature to the best of our knowledge, g = — 1 provides Oi of 0.75 and 02 of 
2.0, independent of the value p assumes, as is obvious from equation- ! 191 since 5 is always multiplied by (g -f 1), which in this 
case vanishes. It is interesting that these a s correspond to p > 2 scaling relations if applied to a p of 2. 

For GRB afterglows, it is not often very easy to infer the value of p unambiguously. The spectral index estimated from 
observations in the optical bands is a composite of the unknown host galaxy extinction and the intrinsic spectral index, 5. 
The X-ray spectrum is not affected by dust extinction but is modified by photoelectric absorption at lower energies. This 
makes the x-ray spectral index to be a function of the unknown gas column density along the line of sight. Also, due to the 
low count rate, it is often difficult to bin the spectrum and get the value of 5 accurately. A third method is to measure the 
flux decay index past the jet break in optical and in x-ray wavelengths and assume it to be p, as predicted by the standard 
afterglow model. Though it suffers from complexities in the modelling of the fireball dynamics, this method is largely followed 
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and trusted. However, the spectral index derived should be consistent with the closure relations between the temporal decay 
index, a and the spectral index, S in various bands. 

Recently several studies have suggested the possibility that the electron energ y index, in f erred by some of the above 
methods, fal ls below 2. Out of the 16 well observed pie-Swift afterglows studied bv lZeh et al.l 1 20061 ) . a2 of five afterglows 
fall below 2. IShen et al.l 1 20061 ) along with blazars and PWNs, study a sample of well monitored X-ray after g lows observed 
by BeppoSAX and Swift. The inferred values of p fall below 2 for eight of them (See figure 5 of IShen et al.l (12006 )). Earl y 
evolution of several x-ray afterglows monitored by Swift have shown an unprecedented 'flat' evolution ( Nousek et al.l 120061 1. 
Though not all of them may have an intrinsic flat electron energy spectrum (some could show shallow decay due to prolonged 
energy injection from the central engine), some are well within the expectations of hard spectrum models. In some of the Swift 
x-ray lightcurves (for example, GRB050820, GRB051109A, GRB061024), the nor mal decay phase, which follows the shallow 
phase, has a values expected from an underlying hard electron energy spectrum ( Liang et al.ll2008l ). 

In the following section, we model three pre-Swift afterglows, with rich multiband data set, showing evidence of an 
underlying hard electron energy spectrum. We consider g as a flt parameter and use a range of —2. < q < +2 while searching 
for the best fit. 



7.1 GRB010222 

GRB 010222 (jPirjboOll ') . at a redshift of 1.477 (| Jha et al.ll2"ooi] : iMirabal et al.ll2002h was one of the first afterglows seen with 



hard electron spectrum and it initiated theoretical work in that direction (BOl, DCOl). 

The optical a f terglow evolution w as initially shallow (ai ~ 0.6) and it steepened to an 02 of 1.3 



Sagar et al 



2001; 



1.4 around ~ 0.5 day 
St anek et al.l 1200 il l. Around the same time the x-ray lightcurve also steepened from ai ~ 0.6 to 02 ~ 1.3 



in't Zand et al.ll200lh . Assuming this early achromatic break to be due to the lateral expansion of the jet, a hard electron 



distribution is required to explain the evolution past thi s break. The spectra l index, 5o, within the optical band was found 
to be 0.89 ± 0.03 after correcting for Galactic extinction ( Mirabal et alJ 2002|'). The x-rav spectral index (Sx) depen ds on the 



assumed value of neutral hydrogen column density of the host galaxy. (|in't Zand et al.ll200ll : lBiornsson et al.ll2002l ). however 
it falls in the range of 0.7 - 0.9. 

Our model with pi ~ 1.5 and q ~ 1.3 reproduces the observed lightcurve decay indices before and after the jet break. 
We assume Vc to be below both optical and x-ray bands at ~ 0.5 day and t^j to be above th e x-ray bands. Along with the 
extinction in the host galaxy {Eb-v ~ 0.03; starburst type extinction law bv ICalzettil (|l997l )) this reproduces the observed 
optical and x-ray spectrum. ^^^^^^^_^^^^^^^ ^_^ 

A model with q of 1.0 and Ui in x-ray bands reproduces the data fairly well ( Bhattacharva fc Resmil 20041') an d also 



explains the spectral steepening seen towards the x-ray band (the x-ray spectrll index derived bv l i^'t Zan TTtn ]^ 



using 

the Beppo-SAX data , is steeper than that in the optical bands). However, our best fit is obtained when q is 1.3, not when it 
is unity. A higher q requires a steeper pi to reproduce the lightcurves decay indices as 5\ and 52 decrease as q increases. The 
best fit with q = 1.3 (figure 5) requires that Vi > Vx- 

We calculated the inverse compton emission for these parameters, and found that it is negligible at the x-ray frequencies. 
We obtain a peak fiux fp of 1.04 mjy and the peak frequency Vm of ~ 200 GHz, at the time of the break. From these fit 
parameters, we infer an isotropic equivalent energy of 5.9 x IO^^tiq erg, a jet opening angle of 2.1°71q , and a total energy 
of 3.6 X IO^^tIq ergs. An upper limit of 10^ is estimated for ^. The best fit model along with the observations are displayed 
in figure 5. The spectral parameters and physical parameters are lis ted in table 4 and table 5 respectively. 

We note that a model assuming continuous energy injection bv lBiornsson et al.l (|2002l ) can also reproduce the observed 
evolution of this afterglo w. Another explana tion for the achromatic break observed around ~ 0.5 day is the non-relativistic 
transition of the fireball ( Masetti et al.ll200lh . but such an early non-relativistic transition would require a very high ambient 
medium density in ~ 10® atom/cc for the observed fiuence of this burst) which would have suppressed the radio flux to 
nano-jansky levels. 



7.2 GRB020813 



GRB020813 was detected by HETE-II (jVillasenor et al.ll2002l l at a redshift of 1.26 (|Price et al.ll2002h . The opt ical aft^ergli 
of th is burst, like GRB010222, exhibited a shallow decay and an early break (ai ~ 0.8 , tb ~ 0.5 day in optical (|Covino et al. 



2003)). The x-ray observations started a fter the optical brea k, the lightc urve exhibited a sin gle power law decay consistent with 
the post break optical decay (ofo ~ 1-4 ( Covino et aL|2003), olx ~ 1.4 ( Butler et al.ll2003l )V The optical photom etric spectral 
index, corrected for Galactic absorption was ~ 0.9 (jCovino et al.ll2003l l and the x-ray spectral index was ~ 1.0 (JButler et al 
20031 1 with no absorption column in excess of the Galactic value of 7.5 x 10^" cm~^. 



The value of p obtained from the best fit model is 1.4, for a g of 1.3. The jet break occurs at around half a day. We 
assumed Vc to be ~ 2.5 x 10^"^ Hz at the time of the break, below the optical bands, to satisfy the observed a and 5 in both 
x-ray and optical frequencies. The synchrotron peak frequency Vm is around 4 x 10^^ Hz at the time of the jet break and the 
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peak flux f^^ is ~ 1.4 mjy. The self absorption frequency Va cannot be constrained using current observations. Our model 
requires additional e xtinction from the host, with rest frame A^ of 0.09 corresponding to an Eb-v of 0.04 and a starburst 
type extinction law (|Calzettilll997l ). 



The derived total energy of the burst is 3.6 x 10 



49„2/5 
^0 



ergs, confined in an opening angle of 2.3°nn . The upper limi t 



1/10 



on 5 is 10 . The polarisation lightcurve of this afterglow has been explained in terms of a structured jet (JLazzati et al.ll2004l l 



The light curve from a stru ctured jet viewed at an angle ^o hardly differs from that of a homogeneous jet with half opening 
angle 9o ( Rossi et al.ll2002l ) (especially for a jet structure described by a 9~^ powerlaw). Hence we can still safely assume the 
shallow powerlaw model for the electron energy distribution within the jet, even though we are not using the structured jet 
calculations. However, The total energy calculations will be affected , if the ener gy distribution is not homogeneous within 
the jet. If we assume that our inferred value of 9o, which accordi ng to Ross^et_alJ will be the viewing angle, is approximately 
equal to the half opening angle of the core of the structured-jet ( Rossi et al.ll2002 'l. and if the actual extent of the jet is 90°, 
the energy inferred will be ~ 9 times smaller than the true energy (see lRossi et al.l for details). 

The best fit model along with the observations are displayed in figure 6. The spectral parameters and physical parameters 
are listed in table 4 and table 5 respectively. 



7.3 GRB041006 

We have presen ted multiband mo delling of this afterglow, which is yet another example of a p < 2 electron distribution, in 
another paper (jMisra et al.ll2005l ). We therefore do not describe this in detail here. We assume the cooling frequency (vc) 
to be below the optical bands to satisfy a of 0.5 and 5 in the range of 0.6 — 0.7 simultaneously. There is no signature of 
steepening seen at the higher energy end of the spectrum from the available observations. Hence we place Vi above the x-ray 
band. We compute the spectral evolution of the afterglow with these basic assumptions. For the sake of completeness, we list 
the spectral and physical parameters from our model in table 4 and table 5. 



8 CONCLUSIONS 

In G RB afterglows, as in other non-thermal sources, the shock accelerated electron spectrum at times assume a hard distribu- 
tion ( Hoshino et al.lll992l : lLeahv et al.lll989r ). But almost all of the theoretical and modelling work in GRB afterglow physics, 
by default, assume a single steep power law for the distribution of electrons in the downstream plasma. The presence of a 
p < 2 spectrum, in a minority of cases, has however not received a fair share of attention. Calculations to derive the physical 
parameters of the burst in such cases are often not done consistently. Early attempts to model GRB afterglows with hard 
electron energy spectrum had several loopholes. 

We have, in this paper, followed the approach of parametrising the temporal evolution of 7^ (thereby leaving room to 
account for different possible physical processes that could determine 7^) as 7^ oc V (BOl) and obtaining the afterglow flux 
decay index for different values of q. We have obtained expressions to calculate the observables from the physical parameters 
of the system which in turn can be used to derive the latter. We present multiband modelling of three afterglows, assuming 
ultra-relativistic expansion, and estimated their physical parameters. 

For all these afterglows, we obtain good fits when g ^ 1. The inferred lower limit of ^ is around 10*. Within the present 
understanding of particle acceleration physics, a mechanism which produces g ^ 1 and ^ ~ 10'* is not known. However, 
future observations of GRB afterglows in the high energy range which can be achieved by upcoming satellites GLAST and 
ASTROSAT will shed more lights on these parameters. For none of the three afterglows, the synchrotron self absorption 
frequency was well constrained. This left us with four observables and five unknowns, so we obtained the physical parameters 
as a function of the assumed value of ambient medium density. Though all of these afterglows were bright in their 7-ray 
output with isotropic equivalent energy in 7-rays ~ 10^^ - 10^^ erg, the total kinetic energy derived from multiband modelling 
is relatively low (~ 10*^ erg). This is partly due to the narrow beaming angle derived from an early jet break (for all the 
jets, 9 is roughly 2.5°). Perhaps kinetic energy being an order of magnitude less than the energy output in radiation could 
be a trait associated with the presence of hard electron energy spectrum. More afterglows and their detailed modelling is 
required to examine this possibility. Another significant characteristic of all the three afterglows is a relatively low value of 
the synchrotron cooling frequency. While for most afterglows discussed in the literature, Vc remain above optical bands longer 
than a day after the burst, the three afterglows discussed here have, in our model, Vc falling below the optical band within 
3 hours. 

The origin of the ha rd electron distribution is not yet clear. Different p hysical processes such as diffusive shock acceleration 
( Achterberg et al.ll200in . cyclotron wave resonance ( Hoshino et al.lll992l ) etc. are beginning to be explored in detail in the 



context of relativistic shocks. Further developments in this area will hold the key to understanding the origin of the observed 
spectra of Gamma Ray Bursts and their afterglows. 
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Table 1. Temporal indices of the spectral parameters. For general q and s 

frequency before jet break after jet break 

3 + (3-6)pi-2q(2-pi)(s-3) 2q-pi-qpi 

™ 2{4-s)(pi-l) 

,, f,, , ,, ^ ,, \ s(10g-4-pi-5pig) + 15(^pi+pi9~29+2) 

l^a{l^a <-ym < Vc) 10(4-s){pi -1) 

,/ (,, ^ „ ^ ,, \ 3(2+pi~4<;+2pi<;)-6pi-20+12,;-6pii3 

ya(,Vm <l'a< Vc) 2(4-s)(pi+4) ^1+4 





Pl-1 






7pi- 


-5pig+10q- 


-12 




10(pi- 


-1) 




(25- 


-4-Pi- 


'9Pi: 


) 



3 + 


22 


5 + 5 
6- 


(s-4) 
-s 


3(s- 
s(l+2q)- 


-4) 
-6(9-H) 


2(4- 


-s) 


3s- 


-4 


2(4- 


-<.) 




s 



Vai^c < l^a < Vm) 

J^i ""^"olA-^r"' -(1 + 9) 



U„ 



2(4-s) 



Table 2. The spectral indices (5) and lightcurvc decay indices (oi; before jet break, Q2; after jet break) for various spectral regimes in 
slow cooling phase. Note that a depends upon the value q assumes. The expressions assume forms similar to those in p > 2 case, if q is 
set to unity. 

spectral segment 5 ai (ISM, WIND) 02 



U < Um < fa < Vc 



(1(1 — 7pi+3pi9— 6q) 6-5pi+pi9— 29 3pi — 6— 3pi9+69 

8(pi-l) ' 4(l-pi) 6(pi-l) 



Ua < 1^ < I'm < Vc 



Um < U < Va 



i P1+P19-29 2-P1+P19-29 -2pi+3-29+9pi 

3 4(Pi-l) ' 6(Pi-l) 3(pi-l) 



5 5 7 

2 4' 4 



Um < U < Ui < Uc 
Vm < U < Uc < Ui 



-^^^ -|(pi+pig-2g),i(2g-pig-2pi-l) _ 2(2^1)^il(l±2) 



vm<ur<u<vc -^^^ -|(p2+P2g-2g),i(2g-p2g-2p2-l) _ HiS^il^liilSl 
Um<i^c<u<u, ~E^ i(6<?-3pi-3pig-2),i(2</-pig-2pi) „ 2(-?-i)-Pi(9+i) 

Vm < I'i < I'c < I' 

-f i(6g-3p2-3p2g-2),i(2<?-2p2-P2g) - 2(9-l)-P2(9+l) 

I'm < Uc < l^i < I' 

Table 3. Same as table 2, but for fast cooling phase. After u goes above both Uc and i/m, the respective positioning of these frequencies 
does not affect lightcurve slope and the indices will be the same as that of the corresponding slow cooling regime. 



spectral segment S ai (ISM,WIND) 02 

i^ Kua <Uc) 2 1,2 1/9 

Ua<u < Uc) 1/3 1/6 , -2/3 -1 

u <Uc<Ua) 2 1,2 13/5 

Uc<u <ua) 5/2 5/4 , 7/4 13/5 

(ua,Uc) <iy <u„^ -1/2 -1/4,-1/4 -1 
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Table 4. Fit parameters of the three modelled afterglows, given around the time of jet break. 



Fit Parameters 


GRB010222 


GRB020813 


GRB041006 


VI 


, 47+0.004 
^•^'-0.003 


1 4n+o-oo7 

-^•^'^-0.004 


1.29-1.32 


P2 


2 n4-°"i 


~2.1 


> 2.2 


q 


1.3 ±0.06 


1.3 ±0.05 


0.95-1.14 


l/m Hz 


2.24+|5-^g X 10" 


3.99+1 j« ^ ,012 


(1.2-3.0) X 1012 


Uc Hz 


9.03o+»3f X 101-^ 


2.33lE!;i^ X 10" 


(1.0-2.0) X lOi* 


Ui Hz 


>10« 


> 5 X IQis 


> 2.4 X 10^0 


fp mjy 


1 n37+o-''i 


, ot.+0.025 
^■•^°-0.065 


(0.37 - 0.49) 


tj day 


n t:c+0.035 
'J-'' ''-0.033 


0.48 ± 0.03 


0.17-0.24 


E{B^V) (host) 


0-035+HgL inag 


0.03t»;0« mag 


0.01 - 0.05 mag 


Host Gal. B band 


25-64l°-^g mag 


- 


- 


" V band 


26.29l°;f mag 


- 


- 


" R band 


25.83l°;f mag 


- 


- 


" I band 


25.59 ±0.25 mag 


- 


- 


" 8.46 GHz 


SSlgA'Jy 


- 


- 


" 4.86 GHz 


SOl^s^Jy 


- 


- 



Table 5. Derived physical parameters for the three afterglows. Since Ua was not well constrained in all the cases, the parameters are 
presented as a function of the ambient density no, normalised to 1 atom/cc. 



physical parameters 


GRB010222 


GRB020813 


GRB041006 


PI 


~ 1.0 


~ 1.0 


~0.8 


<iBnf^ 


n n27+°-°''i 


n. -,+0.004 

^^•1-0.007 


0.07-0.14 


in-^o 


12.0liig^ X 10" 


> 5.7 X lO'' 


> 2.0 X 10"* 


_i 
Eiso riQ ^ ergs 


5.83l?;i^ X 10-« 


3.22t°ill X 1052 


(2.0-4.0) X 1051 


e.ny-To deg. 


2.0° ±0.008 


2.3° ±0.05 


1.7° -2.8° 


2 

Etot rig ^ ergs 


3.60 ± 0.002 X 10*9 


2.21°^ X 10« 


(1.4-3.4) X 10** 
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log({t-to)/day) 

Figure 1. Evolution of spectral breaks Va (top left), Vrn (top right) and Vi (bottom) for different values of q. For comparison, result of a 
single power law with p = 2.2 is also shown (thin line). Vi is not relevant for the p > 2 case however. The parameters used in calculating 
the curves are: ^ = 1, a spherical outflow of isotropic equivalent energy 10^^ ergs and initial lorentz factor 350 in a homogeneous ambient 
medium of density 0.1 atom/cc. The shock microphysics parameters are: ee = 0.1, eg = 0.01, pi = 1.5, p2 = 2.2 and § = 2000. 
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Figure 2. (left top) Sample model optical lightcurve (4 X 10^^ Hz), (right top) x-ray lightcurve (10^^ Hz) and (bottom) radio lightcurve 
for 22 GHz for three different values of q. 
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t > u 




^- t<tj — ► 

Figure 3. The a - S closure relations for various values of q. The left panel shows the closure relations when the observing frequency 
is below the Uc, the right panel is for u > Uc. In the bottom panels a is calculated before jet break. In the top panel, post jet break a 
values are presented. Solid line is for q = 1, dotted line is for q = —0.5 and dashed line is for q = —1. For 5 = 1, the standard p > 2 
scaling is recovered. Note that for q = —0.5, the ininimum possible value of a is 1.75. For q = —1, a does not depend on 5. 
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Figure 4. The predicted compton contribution from hard electron energy spectrum, in comparison with that from a steep spectrum. 
For frequencies less than 10^^ Hz, the contribution from SSA is rather low for p < 2 spectrum. The parameters used for calculation are, 
^iso,52 = 10^1 " = 100, ee = 0.3 and eg = 10"''. For hard spectrum pi = 1.8, p2 = 2.2, q = 1 and ^ = 5000 are used, and for steep 
spectrum a p of 2.2 is used. The displayed spectra are for ~ 5 days post-burst. 
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Figure 5. Multiband model fits for GRB010222. Points : observed data. Solid line : our model, (a) Radio and x-ray lightcurves. The 
4 GHz lightcurve and the 10^* Hz x-ray lightcurvc arc offset by 0.01 and 0.1 mjy respectively for the ease of viewing. The flattening 
seen in radio lightcurves (panel a) are due to the flux of the starburst host SMMJ14522+4301 (see text for details), (b) Optical BVRI 
lightcurves, appropriately offset to avoid clustering, (c) X-ray spectrum at ~ 1 day from BeppoSAX along with the model. 
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Figure 6. GRB020813: Best fit model along with the observations, (i) The top two curves in the left side panel are radio flux in 8.46 GHz 
and 4.86 GHz respectively. For ease of viewing, 4.86 GHz flux is multiplied by 0.01 mJy. The late time flattening in the 8 GHz data is 
not due to the presence of any host. Such flattening is seen in the radio afterglows beyond a few days past the burst, and is suspected 
to be some non-standard behaviour (see Frail et al. 2004) which is not taken care of by our code. The bottom curve in this panel is the 
x-ray lightcurve at 1.2 X 10^** Hz. (ii) The right panel displays multiband optical lightcurves. I band is offset by —5 magnitudes while V 
band is off set by +5 magnitudes. 
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APPENDIX A: CALCULATION OF THE LATERAL VELOCITY OF THE JET 



The a diabatic sound velocity is defined as, Cs = dP/dp where P is the gas pressure and p is the mass density. IChandrasekhai 
1 1939 ) derives the thermal energy density Uk of a mono-atomic gas to be, 



Uk = n 



3K3{e) + Ki{e) 



iKiie) 



2 

mic , 



(Al) 



where n is the particle number density in the gas and mi is mass of a single particle. Q = mic? /ksT, where T is the 
temperature of the gas. Kti{Q) is the modified Bessel function of order n. In terms of temperature, thermal energy density is 
usually expressed as, na{T)kBT, where a{T) parametrises the temperature dependence. It follows from the two expressions 
that, 



a{T) = e 



3K3{e) + Ki{e) 



1 



(A2) 



4^-2(9) 

In the non-relativistic regime, a{T) approaches the familiar value 3/2 and in the relativistic limit, it tends to 3. For a blast 
wave downstream plasma, with single particle rest mass mi, the average thermal energy per particle a{T)kBT can be written 
as (r — l)mic^. i.e.. 



m\c 



3A-3(e) + A'i(e) 



4A'2(e) 



- 1 



= (r-l)mic^ 



(A3) 



from which we identify (3^^3(0) +A'i(0))/4_K'2(6) with F. Temperature of the gas can be solved for, in terms of F by inverting 
this relation. 

But the total energy density is independent of the dynamic regime of the gas and is given by, u = p(? — {Uk + nmiC^)/V 
where p is the total (rest+inertial) mass density. Using this expression we obtain, 

which gives the sound velocity in the downstream in terms of F as, 
r 1^ 1 

i] -^ (^^) 

Let us examine the limiting values of the above expression and check the consistency. In the non-relativistic limit, fcsT <C mic? 
ie., O 2> 1, the Bessel function takes the form 

2 



K^{Q) 



28 



exp (-9) 



1 + 



4n^ 



1 



86 



Substituting eqn. (|A6|l in eqn. (|A5[) : 

^2 _ ksT 4 [l + i|] 

'^ mi 3[H-t] + [l + 4] 

Neglecting terms of the order of 1/0, expression for sound velocity in a non-relativistic gas is reduced to 

2 kaT 

Cs = 

mi 

Now, in the relativistic limit, ie., when 0^1 The limiting expression for Bessel function is, 

1 (n-1)! 



K„{Q) 



2 (©)" 



(A6) 



(A7) 



(A8) 



(A9) 



Substituting the above expression in (12), and neglecting terms 0(Q^), we get for the sound velocity in a relativistic gas. 



2 kgT se 



c 



mi 24 3 (^^°) 

We calculate the the lateral velocity of matter in the fireball as it decelerates, using eqn (|A5|) . When F — » 1, we shift to the 
non-relativistic expression given by eqn 
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